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c3 ■ Abstract 

We express the entropy of a scalar field 4> directly in terms of its space- 
time correlation function W(x,y) = (4>(x)(j)(y)) , assuming that the higher 
correlators are of "Gaussian" form. The resulting formula associates an 
entropy S(R) to any spacetime region R; and when R is globally hyper- 
bolic with Cauchy surface E, S(R) can be interpreted as the entropy of 
(D . the reduced density- matrix belonging to E. One acquires in particular a 

new expression for the entropy of entanglement across an event-horizon. 
Thanks to its spacetime character, this expression makes sense in a causal 
set as well as in a continuum spacetime. 
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As usually conceived of, entropy is local in time, being defined relative to a hypersurface 
E . To compute such an entropy (at least in its "Gibbsian" guise) one must be able to 
identify a density-matrix p(E) which one can plug into the formula S = Trplogp -1 . 
But according to our current understanding, quantum fields are in general too singular 
to admit of meaningful restriction to lower dimensional subsets of spacetime. If it is 
therefore doubtful whether a concept like "state on a hypersurface" can be well defined, it 
is virtually certain that this concept will break down in the discrete context of a causal set. 
For reasons such as these it would be desirable to define entropy in a more global manner 
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by associating it with a spacetime region rather than a submanifold of codimension one. 
A more global notion of entropy would also seem fitting in connection with black holes, 
whose very definition is global in character. 

Of course, black hole entropy will not be understood fully except against the back- 
ground of a theory of quantum gravity and quantum spacetime. Nevertheless it seems 
fair to say that quantum correlations between the hole's interior and its exterior must 
contribute to its entropy. Perhaps, when suitably understood, some such entropy of en- 
tanglement will even turn out to tell the whole story. Be that as it may, a portion of 
the entanglement entropy can be identified without invoking full quantum gravity, namely 
that portion belonging to whatever "matter" fields are present in the neighborhood of the 
horizon. 

When we compute such an entropy in the continuum, we obtain an infinite answer, as 
is well known. Introducing a cutoff £ on the other hand, we obtain S ' en t a ngiement ~ 

A/r, 

which compares favorably with the exact formula S BH = 2nA/£ 2 p , where £ p = V^tcG is 
the rationalized Planck length. [1] We recognize here the familiar area law, but with a 
coefficient that is not easily determined because it is not "universal" . 

On one hand this sensitivity of S to the details of the cutoff is to be welcomed be- 
cause it lets us learn something about the magnitude of £ and the nature of Planck-scale 
physics. (Although they are sometimes viewed as uninteresting, non-universal quantities 
are actually the most interesting if we are after the microscopic physics!) 

But on the other hand, the very fact that the answer is not universal also creates 
a difficulty. When we try to compute S entanglement with finite £ the answer will depend 
on the details of how £ is introduced (sometimes called, rather misleadingly, "scheme 
dependence"). In itself, this ambiguity might not be a problem, but for the fact that it 
tends to call into question the covariance of S. For example in computing entanglement 
entropy for a Schwarzschild black hole, we might think to neglect all contributions from 
within a distance £ of the horizon, but what does this mean? Do we neglect modes within 
a shell Ar = £ where r is the Schwarzschild radial coordinate? Or should we be using 
instead of r, proper distance along the spacelike surface E whose entropy we want (proper 
distance to a lightlike surface being not in itself well-defined)? But then the answer will 
depend on how we extend E away from the horizon H. If for example we define E as a 
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surface of constant Schwarzschild-time t then it will go through the bifurcation 2-sphere. 
But realistic black holes don't have such bifurcation surfaces; and even if they did, our 
attention when we came to consider the crucial question of entropy increase would have 
to focus on portions of the horizon which were farther to the future. 




Figure 1. The black hole horizon H and a hypersurface E whose entropy 
is desired. The dotted line extends E to a Cauchy surface of the full 
spacetime. 

entanglement entropy in a causal set 

The difficulties we have just reviewed make it hard to define an entropy of entanglement 
that is at the same time covariant, but if one replaces spacetime by a causal set the 
difficulties disappear because no cutoff is needed. That, of course, is one reason for thinking 
of the causal set as more fundamental than the continuum. [2] [3] [4] But even if one is 
seeking no more than a "phenomenological cutoff" to render the entropy finite, a causet 
sprinkled into the spacetime in question offers (uniquely as far as I know) a covariant way 
to obtain one. For this reason alone, it is worth asking whether a natural definition of 
entanglement entropy can be found within causal set theory. 

Let us consider then, the problem of defining a horizon-entropy in a fixed, background 
causal set. First of all one must attach a meaning to the words "event horizon" . Although 
this might seem to pose a problem because a causet contains no subset that could play the 
role of a null surface, the distinction between the interior and exterior of a black hole still 
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makes sense because it depends solely on the causal order. Since quantum entanglement 
concerns only the interior and exterior regions, and not the horizon per se, there is in fact 
no problem. 

The next step, if we were to mimic the continuum story perfectly, would be to define 
the concept of a spacelike (or achronal) surface E , after which we would try to associate 
with each suitable E some sort of "momentary state" or density matrix p E . In fact, there 
does exist a natural choice for E , since we can define it as an antichain, or better a 
maximal antichain (or "slice"). That which is not readily available in a causet, however, 
is anything like p E : a causet-based field theory admits no obvious notion of "state on a 
hypersurface" . 

More generally, none of the surface-quantities which one routinely defines in the con- 
tinuum seems to carry over naturally to an antichain, except possibly a measure of (spatial) 
volume. In fact, it seems to be a rule of thumb that "purely spatial" concepts do not live 
happily in a causet * including for example the concepts of induced metric g^ k and extrin- 
sic curvature Kj k ~ gj k . For this reason a causet is not congenial to canonical quantum 
gravity, albeit one can define many surface-quantities in an approximate sense, including 
both metrical ones and topological ones like homology groups [5] . 

The question then is whether we can free S entanglement from reference to a density 
matrix localized to a hypersurface. In seeking the answer, I will limit myself to what is 
perhaps the simplest case, that of a free scalar field 4>{x) . For this case a corresponding 
theory exists in the causet, and it has been given both algebraic and path-integral (or 
"quantum measure" ) formulations. Inasmuch as the path-integral formulation extends to 
the interacting case, it is arguably preferable to the operator formulation. Nevertheless, I 
will work in the latter context, first of all because more tools are available there, and more 

t More generally, there is evidence for the rule that it is the continuum concepts which 
are "C°-stable" that have simple causet counterparts. 
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importantly because we seem to lack any definition of entropy couched in the language of 
histories. b 

Could it be that this lack betokens an inherent dependence of entropy on some notion 
of "state localized in time"? Were that true, the search for a histories-based definition 
of entropy would necessarily be futile, but such pessimism is called into question by the 
global definition of entropy we will arrive at below in the algebraic context. If indeed 
a more (Lorentzian and non-equilibrium) histories-based concept of entropy is out there 
somewhere, its discovery would, to my mind, mark important progress in the path-integral 
formulation of quantum theories. 

Entropy from the spacetime correlator 

Let us recall some possible definitions of entropy in the context of quantum field theory. 
Conceived algebraically (as it commonly is) , a quantum field is a collection of operators * 
4>{x) acting irreducibly in a hilbert space ij , together with a "global state" represented by 
another operator p (called "the density-matrix" ) such that Tr pA yields the "expectation 
value" (A) of A . One can then define the entropy as S 1 = S(p) = Trplogp -1 . (In 
practice the p that enters this definition will seldom be the "exact microscopic density 
matrix". Rather it will be some coarse-grained version thereof whose entropy is nonzero 
even when the underlying "microscopic state" is pure.) 

Now let R be an arbitrary spacetime region. The operators <f>(x) for x E R generate 
a subalgebra 2l fi of 21, and the expectation-value for 21 restricts to a similar functional on 

b I'm ignoring here the trick that computes the thermodynamic partition function via a 
Wick-rotated path integral. That procedure is indeed histories-based, but the histories in 
question develop in imaginary time, and more importantly, the entropy defined thereby is 
limited to states of thermal equilibrium. 

* In the continuum the symbol <f>(x) must of course be interpreted formally since the 
field 4> is at best an operator-valued distribution. In the following lines I will ignore 
this technicality and whatever difficulties flow from it. They are peculiar to continuous 
spacetime and do not seem to be germane to the causal set context of primary interest in 
this paper, within which the symbol <p 3 belonging to element j of the causet will literally 
denote an operator, albeit an unbounded one. 
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2l fl . In general, the subalgebra 01 R will no longer act irreducibly in ij. Assuming however, 
that we can represent %l R in some other Hilbert space S) R in which it does act irreducibly, 
and assuming further that we can find in Sj R a density matrix p R such that (A ) = Tr p R A 
for operators A e 21^, we can go on to define S(R) = — Tr p R \ogp R , the entropy of p 
relative to the region R. 

Remark Let 21 be the algebra generated by the 4>{x). As such, it is a concrete algebra 
of linear operators in f), but in what might be called the "strictly algebraic" formulations 
of quantum field theory, 21 is regarded as an abstract ★-algebra, while a density-matrix 
p is regarded simply as an expectation-value-functional on 21. Since the set of all such 
functionals is convex, any given p will be (either exactly or to a close approximation) a 
convex combination YlaPaPa °f extremal (or "pure") states p a . In this abstract set- 
ting, S(p) could perhaps be defined as the infimum^ over all such sums of the quantity 
SaPa logfa 1 • I n this way, one could imagine defining the entropy S(R) directly in terms 
of the expectation-value functional, without recourse to either S) R or p R . 

The definition we have just given of the entropy of an arbitrary spacetime region R 
will be the basis of the causal set construction we are seeking. First though, we need to 
relate the entropy of a region to the entropy of a hypersurface. To define the latter itself, 
we can imagine repeating the same steps that led to S(R), only with R shrunken down to 
a hypersurface E , and with 2l E then being the algebra generated (formally) by the "initial 
data" operators <j)(x) and <j)(x) for x G E . 

The point now is that when E is a Cauchy surface for the region R , we have 2l E = 2l H , 
and therefore S'(E) will coincide with the more globally defined entropy S(R). By taking 
for R the so-called domain of dependence of E , we can thus express the entropy of any 
desired hypersurface as the entropy of a spacetime region. And this in turn will let us 
express it directly in terms of the correlation function (<p(x)(/)(y) ) . 

In application to a black hole spacetime, E would be the exterior portion of the 
hypersurface for which we desired the entanglement entropy, and the region R would then 
be given by 

R = future(E) n past (if ) = £> + (E) (1) 
t For dimfj < oo this is proven in [6]. 
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Remark Time-reversing the above, we could also have taken R = D~(E), but the 
analog of (1) would not hold with that choice, making it apparently less convenient when 
transposed to a causet. 




Figure 2. The entropy of the surface £ can be identified with the entropy 
of the region R 

In summary, entanglement entropy is a special case of the entropy S(R) of a spacetime 
region. Our task now is to find a simple formula for S(R) , assuming, as always, that is 
free and "gaussian". 

Sketch of a derivation of a formula 

What follows is in parts still a "work in progress" ; not every loose end has been tied up. 

We have arrived at an algebra 21 = 21(1?) and an expectation-functional or "state" 
(■ ) thereon. In order to derive from these structures an entropy S(R), we will seek to 
represent 21(1?) irreducibly in a hilbert space fj = S)(R) and to find therein a density- 
matrix p : — > Sj such that (A ) = Tr pA for all A e 21. In what follows I will omit the 
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reference to R, as if R were the full spacetime or causet as the case may be. There is no 
loss of generality in doing so, since we will not need to refer again to spacetime points 
(respectively causet elements) outside of R. Also I will write u (p Jn for an individual field 
value. This notation is more convenient than u (p(x) n and it is also more apt for a causet 
C , where the index j runs literally from 1 to N , N = \C\ being the cardinality of C. 

Remark The condition that 21 act irreducibly in S) is crucial here. If we dropped it, 
we could, by "purification", always find a representation in which Trplogp would vanish 
because p would be the projector onto a single vector in $). This would be the case in the 
so-called GNS representation, and for that matter, in the representation of the 4> 3 in the 
original, global hilbert space, assuming the original global state to have been pure. In view 
of these observations, it's somewhat disconcerting that, strictly speaking, the existence of 
an irreducible representation is not always guaranteed. Consider for example the "region" 
of a causet consisting of a single element e. The corresponding subalgebra 21 will be 
generated by a single field-value 0(e) =: q. Now by assumption, (q) = , but in general 
(qq ) will not vanish. If it does not, then q cannot act irreducibly in any representation 
(because then Schur's lemma =>- q = c G K. =>- q = (q) = =>- (qq) = (0) = 0). In such 
a case, the entropy could not be defined by the means we have adopted. However, the 
"strictly algebraic" definition we remarked on earlier would still seem to apply and to yield 
5 = 00. In the following, we will simply assume that an irreducible representation does 
exist. The formula (9) which we will derive will not actually depend on this assumption 
if interpreted as in (8). It would however yield S = in our one-element example, rather 
than S = oo. 

Because we are dealing with a free field, the commutator of ft with 4> k is a c-number 
and we can write 




= iA j 



(2) 



for some real, skew matrix A jfc . 



We define similarly the Wightman function as the hermi- 



tian matrix 



W jk = (<fr> > fc ) . 



Notice that A is twice the imaginary part of W . 
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Assume now that our theory is "Gaussian" (or "Wickian") in the sense that the 
bosonic form of Wick's rule holds with {</>) =0: 

(#-..0) = J](00)...(#) 

(The sum is taken over all ways of pairing the fts, and within each pairing the original 
odering must be preserved.) Thanks to this assumption the "two-point function" W jk 
determines the theory fully. 

For example its imaginary part determines the "equations of motion" for <p. Thus 
(thinking for a moment of the theory on the entire spacetime or causet), if the sequence 
ctj is in the kernel of A in the sense that = then A = ^ctjft vanishes as 

well, (proof: When a G ker A, A commutes with every ft, as follows from the calculation 
[ft, A] = [ft ,a k (f) ] = iA° ' a k = . But since the ft act irreducibly in ft, this implies 
that A = cl, and then c = (A ) = otj{ft ) = whence A = .) This last fact yields a set 
of linear relations among the operators <fi. In the continuum they are just the equations 
of motion for <f>(x) (in our case the Klein-Gordon equation), but in the causet they yield 
only relatively few conditions because the exact kernel of A is typically rather small. 

Now our algebra 21 is generated by the individual field- variables ft , and we can there- 
fore characterize it fully by specifying the relations among these generators. But since 
we have taken our field to be non-interacting, this is quite simple to do. In addition to 
the linear relations just described there are only the bilinear relations (2) (the "canonical 
commutation relations"). Moreover, we have just seen that any zero-eigenvector of the 
commutator-form A lJ results in a linear dependence among the ft . Hence, by passing to 
a linearly independent set of generators, we can assume that A lJ is invertible. 

It follows that there exists a basis for the ft of the form q a , p a (a = 1 • • • n) such that 

By construction this basis block-diagonalizes A, meaning it block-diagonalizes the skew 
part of W . (Recall that W jh - W kj = iA jh = [ft ,4> k ] .) I claim one can also choose it to 
diagonalize the symmetric part, 

R jk = ReW' k = ^({ft,cf> k }) , 
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hence to block-diagonalize W^ k itself. b With this, our problem splits up into a product 
of individual problems, each concerning a single pair (q, p) — a single degree of freedom 
- that you might imagine as position and momentum for a free particle or harmonic 
oscillator (or as the field-components (p J for a 2-element causet, a 2-chain). It thus suffices 
to evaluate the entropy in that special case. 

the entropy for a single degree of freedom 

Our task is now as follows. Given a conjugate pair of variables q and p such that [q,p] = i , 
and given the correlators 

(qq ) (pp ) Re (qp ) 

for a Gaussian density-matrix p , to find S(p) = Trplogp -1 . That p is Gaussian means 
that in a (/-basis it takes the form 

p(q, q') = (q\p\q> ) = (est) exp (-|(, 2 + ,' 2 ) + ^ (, 2 - ,' 2 ) - ^(q - q'f) 0) 

Then S(p) must be some function of the parameters A, B, C, and our task is to determine 
this function. 

To that end, notice first that S must be dimensionless and invariant under unitary 
transformations. From this it can be shown that S can depend on the correlators only in 
the combination (qq ) (pp ) — (Re (qp )) 2 = det Rj det A , where in this simple 2x2 situation, 
A and R reduce to 



A-2Ima e ^^ ) {qp) ] - ( ° 1 
A - 2Imag {(pq) (pp))-{-l 



and 

f (qq) (qp) \ = ( (qq ) Re (qp ) 
V (PQ ) (PP) ) V Re (QP ) (PP ) 



R = Re 



b In fact, diagonalization is possible in general, without assuming that the <\P are linearly 
independent. Given any positive-semidefinite matrix W ', with associated real matrices R 
and A such that W = R + iA/2, one can introduce a (real) basis in which R is diagonal 
and A is block-diagonal with 2x2 blocks. Moreover one can arrange that the non-zero 
matrix elements R 3J = 1. The proof is slightly too long to reproduce here, but see [7] for 
a similar result that applies when A is invertible. 
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It turns out that (as a short calculation will confirm) 



(qq)(pp) - (Re(qp)) 2 = C/2A + 1/4 

whence the entropy depends only on the ratio C/A, while B drops out entirely. But with 
B set to 0, we can take over from [8] the entropy calculation there, with the result * 

_ //log// + (l -//) log(l-//) 
1 — /i 

where 

_ y/1 + 2C/A - 1 

^~ y/1 + 2C/A + 1 

Putting the pieces together yields now 

S = (a + 1/2) log(<r + 1/2) - (a - 1/2) log(a - 1/2) (4) 

where I have defined the eigenvalues of A -1 /? (which are purely imaginary) to be ±ia : 
spectrum(A _1 i?) = ±ia . 

This nicely symmetrical expression is already quite simple, but we can simplify it still 
further by swapping the eigenvalues of A _1 i? for those of A _1 VF = A -1 /? + i/2 . (Recall 
that W = R + zA/2.) Writing these latter eigenvalues as ±iu ± = i(| ± a) , brings S 
finally to the convenient form 

S = u + loga;_|_ — u_ logo;_ (5) 

Notice in this connection that positivity of the correlation matrix (also known as the 
uncertainty principle) implies that a > 1/2 , whence u + and u_ = u + — 1 are both > . 

the entropy in full 

The formula (5) belongs to a single degree of freedom, corresponding to one block of our 
block-diagonalized matrix W. Summing (5) over all the blocks then yields the total entropy 

* Might there be a more direct route from (3) to S(p), perhaps via a diagonalization of p 
or via a use of the "replica trick" ? 
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in the form of a sum over the spectrum of the full matrix A W , which for short I will 
call L, or more conveniently iL in oder that its eigenvalues be real: 

A~ 1 W = iL (6) 

Denoting the eigenvalues of L by A, we have then 

S = ]TAlog|A| (7) 

We have seen that the eigenvalues A are all real, even though L is in general neither real 
nor symmetric. We have also seen that each negative A is paired with a positive eigenvalue 
1-A. 

As derived, this last formula relies on the fact that the q a and p a are all linearly inde- 
pendent. When we revert to the original matrices W-* and A jfc , we will encounter (after 
their diagonalization) a block submatrix consisting entirely of zeroes. In this subspace, 
A is obviously not invertible, but neither is there any contribution to S. To adapt our 
formula (7) to this circumstance, we can rewrite the eigenvalue equation for A _1 VF so that 
it ignores this block of zeroes. To wit, we define the eigenvalues A as the solutions of the 
equation W^ k v k = iXA^ v k , where it is understood that A-* v k must not vanish (just as v 
must not vanish in the more traditional eigenvalue equation Lv = Xv). In other words, we 
define the eigenvalues of L as the solutions of the equation 

Wv = iXAv (Av ^ 0) . (8) 

With this understanding (and with the usual convention that OlogO = 0), we have finally 

S = Tr L log \L\ = ^ A log |A| (9) 

A 

A remarkably simple formula! 

Remark Derived under the assumption that the cjp can be represented irreducibly in 
some Hilbert space, equations (8) and (9) give a good account of the entropy when the 
kernels of A and R (or equivalently W) coincide. From the condition W > that W 
is positive semidefinite, it follows that kerA D ker R, but the converse inclusion is not 
guaranteed in general. When it fails, no irreducible representation exists (as seen in an 
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earlier remark), and our entropy is to that extent ill-defined. If we adopt the more general 
"strictly algebraic" definition described earlier, the resulting entropy will diverge. By 
way of illustration, consider our earlier example of a single "field operator" (p 1 = q with 
(q) = and (qq) = s 2 (and with 21 being the algebra of polynomials in q). One can 
think of q in this case as a classical random variable with gaussian distribution function 
1/V27rs 2 exp(— x 2 /2s 2 ), the pure states corresponding to specific sharp values q = x for 
x G R . Since these states are uncountably infinite in number, our "strictly algebraic" 
entropy will be infinite, absent some sort of short-distance cutoff or regulator. 

Remark Given the commutator function A jfc , the construction of [9] and [10] produces a 
distinguished "vacuum" by taking the two-point function W to be the positive projection 
of iA, i.e. by setting W = ^2' n vj n \n)(n\, where the sum is over the postive eigenvalues 
w n of the matrix iA = J2 n w n\ n )( n \- Evidently, the positive eigenvalues produce solutions 
of (8) with A = 1, while the negative eigenvalues produce solutions with A = 0. The net 
entropy therefore vanishes, as it must since the vacuum in question is the "ground state" 
of a Fock-type representation, hence a pure state. 

Research at Perimeter Institute for Theoretical Physics is supported in part by the Gov- 
ernment of Canada through NSERC and by the Province of Ontario through MRI. 
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